Self similarity and attraction in stochastic 
nonlinear reaction-diffusion systems 



Wei Wang* A. J. Roberts^ 
6 November 2011 



Abstract 

Similarity solutions play an important role in many fields of sci- 
ence: we consider here similarity in stochastic dynamics. Important 
issues are not only the existence of stochastic similarity, but also 
whether a similarity solution is dynamically attractive, and if it is, 
to what particular solution does the system evolve. By recasting a 
class of stochastic PDEs in a form to which stochastic centre mani- 
fold theory may be applied we resolve these issues in this class. For 
definiteness, a first example of self-similarity of the Burgers' equation 
driven by some stochastic forced is studied. Under suitable assump- 
tions, a stationary solution is constructed which yields the existence of 
a stochastic self-similar solution for the stochastic Burgers' equation. 
Furthermore, the asymptotic convergence to the self-similar solution is 
proved. Second, in more general stochastic reaction-diffusion systems 
stochastic centre manifold theory provides a framework to construct 
the similarity solution, confirm its relevance, and determines the cor- 
rect solution for any compact initial condition. Third, we argue that 
dynamically moving the spatial origin and dynamically stretching time 
improves the description of the stochastic similarity. Lastly, an appli- 
cation to an extremely simple model of turbulent mixing shows how 
anomalous fluctuations may arise in eddy diffusivities. The techniques 
and results we discuss should be applicable to a wide range of stochas- 
tic similarity problems. 
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1 Introduction 

Consider the stochastic field u{t, x) governed by the nonlinear react ion- diffusion 
stochastic partial differential equation (spde) 

lit = n^x + f{d^, u, + gid,, u, t-^'^)B{t, x) (1) 

on an infinite spatial domain in one dimension. Here the drift nonlinear- 
ity /((9x, u, is cubic in the field u, the spatial derivative, and decaying 

time stationary examples include the cubic reaction / = — u^, the 



self-advection nxix of Burgers' pde, and the nonlinear diffusion esign{xixx)iixx 
explored by Barenblatt [3^ §3.2.1]. For some stochastic process B, with char- 
acteristics defined variously in later sections, the noise term has nonlinear 
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coefficient g{dx,u,t~^/'^) which is quadratic in the field u, the spatial deriva- 
tive, and for example, the stationary linear advection u^- Equation ([T]) 
informally describes the scope of spdes we consider. 

Interpret all stochastic calculus in the Stratonovich sense so that, with 
care, the normal rules of calculus apply. 

We apply some stochastic theory together with centre manifold theory to 
help understand and solve the long time evolution of such stochastic diffusion 
with nonlinear reaction-advection. But on the infinite spatial domain there 
is no clear cut slow eigenspace for spde ([T]). For example, in the absence of 
noise and nonlinearity, f = g = , one substitutes spatial structures with 
wavenumber k into the diffusion pde ([T|, u oc e*'^^'+^*, to find the continuous 
spectrum A = —k"^. Such continuous spectra do not have a well defined 
decomposition between fast transients and the slow, long lasting, modes of 
interest. 

A special transformation changes our view of the dynamics of spde (jT| 
into one with a clear fast-slow separation. We extend to stochastic dynamics 
the transformation used so interestingly by Wayne et al. [U HH HT] for deter- 
ministic similarity, and analogous deterministic analysis on appropriate ini- 
tial conditions Introducing log-time and similarity variables transforms 
the stochastic problem to one of seeking a stochastic field m(t, ^) where 

X 1 
r = logt, ^=^' ^ " 7!^*^^'^''' 

Then the dependence upon the scaled space variable ^ causes the diffusive 
Gaussian spread following a point release. 



a 

exp 



2Vvrt 



(3) 



to correspond to a fixed point of the dynamics for the stretched field u, 
namely 



a 

u* = ^z—^ exp 
Z\/tt 



(4) 

Further, we shall see that the algebraic decay in real time t from any compact 
release to the Gaussian ([3| transforms to an exponentially quick decay in log- 
time r to the fixed point Q. Solutions of the spde ([T]) that are approximately 
the self-similar Gaussian spread will be approximated based upon analysis 
about the fixed points ^ in log-time 

A similar, but stochastic, version of the transformation ^ shows the 
appearance of an eddy diffusion model of the long time behaviour of a sim- 
ple turbulent mixing model introduced by Majda, McLaughlin, Camassa et 
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al. [321 |33l [m [13] . Whereas they primarily explored solutions statistically 
stationary in space, as is the theme in this article, Section |5] characterises 
the spread following a compact release. Although a mean eddy diffusion 
emerges, significant anomalous diffusion occurs throughout due to stochastic 
fluctuations. 

Section |2] proves the existence and emergence of stochastic self-similarity 
for a stochastic Burgers' equation in the class of the spde ([T]). Burgers' 
deterministic partial differential equation for a field u{t, x) is the case / = 
—uUx and g = oi the spde ([T|, namely Uf = Ua:x — ■ This pde was 
proposed by Burgers' [I2] to illustrate the statistical theory of turbulent fluid 
motion. To better model the turbulence and turbulent flow in the presence 
of random forces, stochastic Burgers' equations have been suggested [HI [T5|, 
123 [2ni EHl iSl SS] and studied by many people [SlElIITllSniESlEHlEniES]. 
We explore the stochastic solutions u{t, x, u) to the stochastic Burgers' pde 

Ut = u^x - uu^ + x) (5) 

where ^{t,x) is some stochastic force, to be detailed later, defined on a 
complete probability space P). Figure [l] shows just one realisation of 

solutions to Burgers' spde ^ to illustrate the emergence of stochastic self- 
similarity. On bounded domains Da Prato et al. [38] proved the existence 
and uniqueness of global solution when the noise term Q3(t,x) is a noise 
white in time and fixed spatial structure. Holden et al. [23] also derived the 
same result by some white noise calculus. On an unbounded domain Bertini 
et al. [5j constructed a global solution by a Cole-Hopf transformation with 
space-time white noise x). 

Here we consider a family of solutions with special spatio-temporal form, 
stochastically self-similar solutions, of the stochastic Burgers' equation ([s]) on 
the unbounded real line with a particular form for the stochastic force 23 (t, x). 
The existence of self-similar solutions and the asymptotic emergence of a 
self-similar solutions. Figure [T| describes the self-similarity of the stochastic 
Burgers' equation. 

Self-similarity is an important property of some convective diffusion equa- 
tions, of which Burgers' equation is a special case. Many researchers have 
established the existence of self-similar solutions of deterministic systems [TOl 
|22| [311 [m [50] , and described the asymptotic behaviour of such self-similar 
solutions [H [22], |35l [50] . Here for stochastic partial differential equation, we 
construct a self-similar solution in the sense of distribution, which is called 
a random self-similar solution. Previously, Eyink and Xin [2T] studied the 
deterministic similarity of statistical quantities. In contrast, here we addi- 
tionally explore the structures of the stochastic fluctuations. 
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Figure 1: movie of stochastic self-similarity emerging in a realisation of a 
stochastic Burgers' equation. Here simply obtained via a naive but fine scale 
discretisation of the spde (|5| for noise ?B oc W{t,x)/t. The top frame show 
the decay and spread in physical variables, whereas the bottom frame shows 
the same realisation in similarity variables ^ and hence the approach to a 
stochastic self-similarity. 



5 



1.1 A stochastic slow manifold emerges 

Let's look at the transformation ^ applied to the general react ion- diffusion 
SPDE ([!]). First, consider the coefficient functions: because the reaction 
term / is assumed cubic in its arguments f{dx,u,t~^^'^) = t^^/^/(9^, m, f ); 
similarly, as g is assumed quadratic in its arguments the noise coefficient 
g{dx,u,t~^^'^) = t~^g{d^,u, 1). Second, the original SPDE Q thus transforms 
to the SPDE 

Ur = Cu + fid^, u, 1) + g{d^, M, 1)W, (6) 

where the new noise W = y/tB is assumed to be Q- Wiener space-time noise, 
and where the linear operator 

Cu = u^^ + ^^u^ + . (7) 

The long time dynamics of the original spde Q are revealed by the dynamics 
of the transformed SPDE ([6]). 

Consider the SPDE (|6| for negligible noise, W = . The SPDE then has 
equilibrium u = . Linearised about this equilibrium the SPDE has discrete 
spectrum = —k/2 , k = 0,1,2, ... , with the corresponding eigenfunctions 



efc(0 = CkHk{,i/V2)G{C), for Gaussian G{i) = — ^exp 



-e 
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(8) 



where the Hermite polynomials -f^fc(C) = (— l)'^e^^/^(9'^e~'»^''^/9C'^. These 
eigenfunctions forms a standard orthonormal basis: we choose the normalisa- 
tion constant Ck = {2y/7r/k\Y^'^ as then ek{^Ye^ = 1 . The wonder- 
ful aspect of this transformation to log-time is that the spectrum of the linear 
operator £ is discrete and there exists a clear separation of the fast modes, 
mode numbers k > 1 , from the slow mode, k = . That is, the equilibrium 
u = of the SPDE ^ has the slow subspace u = aG{^) for all a. We use the 
existence of this slow subspace as the basis for analysing emergent stochastic 
self-similarity. 

We suppose the noise is Q- Wiener with the following spectral decompo- 
sition 

oo oo 

^(^'0 = J2^kWk{r)ek{0 = J2^kMr)Hk{^/V2)G{0, (9) 

k=0 k=0 

where Wkir) are independent Wiener processes, the noise coefficients hk decay 
sufficiently rapidly with mode number k (for previous rigorous theory only 
a finite number of noise coefficients hk can be non-zero), and the overall 
magnitude of the noise terms are denoted by some norm b = \\bk\\. Because of 
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the spectrum, the conditions on the noise, and for smooth enough coefficients 
/ and g, there must exist a stochastic slow manifold [H [1]. 

Furthermore, under suitable conditions on the noise and the SPDE coeffi- 
cients / and g, the stochastic slow manifold is exponentially quickly attractive 
to all nearby initial conditions. Thus stochastic self-similarity of solutions to 
the SPDE ([T]) emerges from generic compact initial conditions. 



2 A stochastic Burgers' equation 



This ffist extensive section proves in detail the emergence and nature of 
stochastic self-similarity in the stochastic Burgers' equation (|5|. Subsequent 
sections on other spdes are less rigorous but apply more widely. 

We here construct and analyse a stochastic self-similar solution for the 
stochastic Burgers' equation (|5| for t > 1 with a special stochastic force rj. 
For this we invoke the self-similarity transform (|2]), then 



du = [ugg + ^^u^ + — uu^] dr + dW{T,^) 



(10) 



with initial u(0,^) = u^{^) = u{l,C,). Stationary solutions u{^,u) to equa- 
tion (10) is a self-similar solution of stochastic Burgers' equation ([s]). Here in 
order to construct a self-similar solution of stochastic Burgers' equation ([s]) 
we assume that stochastic process W{t,^) is an L^(M) valued Q- Wiener pro- 
cess defined on {Q, J-", P) with covariance operator Q which is detailed later. 



To construct a stationary solution of (10), we consider the system in a 



weighted space L'^{K) defined in the next subsection 2.1 First by using en- 
ergy estimates and the compact embedding results of the weighted space, we 
show the tightness of solution with initial value in the space L'^{K). Then 
the classical Bogolyubov-Krylov method pQ implies the existence of station- 
ary solution of (10). Further, we show the locally attractive property of 



the self-similar solution by a Cole-Hopf transformation (Theorem 13) and 
a local random invariant manifold method (Theorem 18). The Cole-Hopf 



transformation makes the equation ( 10 ) to be a linear one which shows that 



the stationary solution is determined uniquely by the mass of the solution 
provide the stationary solution is small enough in the space L'^{Q, L°°{M.)). 
Then every small solution in the space L^(f2, L°^(]R)) is attracted by a unique 
stationary solution. This is also shown by a local random invariant mani- 
fold discussion without assuming solution is small in L°°(M). The last sub- 



section 2.5 shows the globally asymptotic convergence in probability which 



shows the existence of global random invariant manifold of equation (10). 
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2.1 Preliminary 



Consider the stochastic pde (10). For this we use the hnear operator L 
defined in equation ([T]). Define the exponentially increasing weight function 
-^(0 — 6xp{^^/4}, and then introduce the following weighted functional 
spaces for j9 > 

and for positive integer k 

H\K) = LeL'iK):\\u\\l,^^^= WuWIhk) < oo 

I 0<a<fc 

Then linear operator C is self-adjoint and generates an analytic semigroup S{t) 
on the space L'^iK) with the domain D{C) = H'^{K) gH]- Recall that the 
eigenvalues of operator C are = — A;/2, k = 0,1,2,..., with the corre- 
sponding eigenfunctions ^ forming a standard ortho normal basis of L'^{K) 
as ||efc||L2(;^) = 1 . 

Define inner product on space L'^{K) as 

{u,v)= [ uiOviOKiOd^, u,vEL\K). 
Jr 

Then denote by Pq and Pg the linear projection from L'^{K) to the slow 
subspace Eq and the stable subspace Es respectively. Then Eq = span{eo(0} 
and 

E, = E^ = G L\K) : J^^iO = o} • 

This first lemma presents some basic properties on these weighted spaces 
Lemma 1. 1. The embedding H^{K) C L'^{K) is compact. 

2. There exists C > such that for any u G {K) 

[ \uiO\'\eK{Od^<C [ iVuiOl'KiOd^. 
Jr Jr 

3. For any u G H^{K), 

\um'K{Od^< I \Vum'K{i)di. 
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4- For any u G Eg , 

{Cu,u) < • 

5. Ifue H\K), then K^/^u G 

6. For any q > 2 , e > there exists constants C^m > , i? > , such that 



for any u G H^{K) fl L 



q 

loc 



||m|Il2(^) < e||^?||L2(x) + C'e,<jll^lli9{B{0,R)) 



Remark 2. By part |3] in the above lemma, in the space H^{K) we can define 
norm ||Vm||l2(x) which is equivalent to . 

Further, by the spectrum property of the linear operator C we define (£ — 
1/2)'^ for any 7 G M [49J. Then define Sobolev space H'^{K), for any 7 G 
R as P((£ - 1/2)^/2), the domain of (£ - 1/2)^/2, and by the embedding 
theorem |19] , H'^^ is compactly embedding into H'^'^ for 71 > 72 . 

For our purpose we assume that under the self-similar variable the special 
stochastic force 77 is assumed to be an L^(A')-valued Wiener process W{t,^), 
r > , which is conservative; that is, it satisfies 



/ 



W{T,Od^ = 0. (11) 

Moreover, W{t,C,) has the following spectral expansion 

00 

W{T,0 = J2^kWk{r)ek{0 (12) 

k=l 

where {wk{T)}k are mutually independent standard scalar Wiener processes. 



By the conservation assumption (11), 60 = and so is omitted from this 



sum. So W{t,^) = Ws{t,^) the part of W{t,^) in the stable subspace Eg. 
Furthermore we assume 



J2^l<oo. (13) 



k=l 



For later purpose we extend W{t,C,) to the whole time interval (—00, +00) 



by W{—T,^) = —W{t,^), r > . In the following we consider equation (10) 
on the canonical probability space (^o, J-o,IPo) which consists of the sample 
path of W{-,u) in the space C{R, L^{K)), and denote by {9^-}^ the Wiener 
shift on [Ij. 

Remark 3. The special assumption on B(t,x) does not exclude the existence 
of self-similar solution for other cases. 
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Remark 4. Notice here that fio is different from fi. However by the self 
similar transformation, for any u & Q, there is a sample path of W{t,C,), 
then there is a cjq G Qq represents one sample path of W{t,^) in Qq . The 
inverse is same. So in the following we see Qo same as Q . 



To estimate the solution to (10) we need the Ornstein-Uhlenbeck pro- 



cess rj°'{T,^) which solves the following equation 

rfr/" = [£r/" - ar^"] dr + dW{T, 0, (14) 

for some a > ; that is, f]"{T,^) = e^^~°'^^'^~'^^dW{s). (Here we introduce 
the differential dW{s) to denote the log-time differential dW{s,^) but with- 
out explicitly including space ^ to denote the differential does not involve ^.) 
Then we have the following estimates on ?7"(t, ^). 
Lemma 5. For any e > and p >2 there exist a > such that 



Proof. We just need to prove this result for p = 2j , j > 1 



First the unique stationary solution to equation ( 14 ) is written as 



r/"(r) 



5^ A / e(-^^-")(^-^)ciu;fc(s)efc 
k=i •^-'^ 



Then by the property of stochastic integral, 

1 2 \ _ bkXk 



E(hFMIlW))=E 

k=l 



By the assumption (13) on iy(r, ^), Lemma [T] and Remark |2] we have the 
result for j = 1 by choosing a large enough. For j > 1 , the result is followed 
by the Guassian property of rj" [STJ Lemma 7.2] . □ 



Remark 6. Here the a is chosen to be large so that r^" is small. Subsection 3.1 
shows, for a = , that the process r/" is the emergent stochastic slow manifold 
of linear system ([l4]) with a = and describes the stochastic self similarity 
of the linear problem ([T]) with / = and g = 1/t . 



Let T]" = t]q + r]'^ G Eq (B Es , then by the assumption (11 ), ?7q = . By 
Lemma [T] and [5| we then have 

Corollary 7. For any e > 0, there exists a > such that for any p > 1 
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2.2 Existence of a stochastic self-similar solution 



For any r > , in the mild sense equation (10) is written as 



u{t) = S{t)u^+ S{t- s)u{s)ui.{s)ds+ S{T-s)dW{s). (15) 
Jo Jo 

Then by the theory for abstract stochastic evolutionary equations [37] we 
have the following theorem. 

Theorem 8. For any T > and m° E L'^{no,L^{K) n L°°(]R)), there is a 



unique solution u{t,C,) to SPDE (10) in the space L'^{Qq,C(0,T; L'^{K)) fl 
L\0,T;H\K))). 

Proof. Notice that the nonlinearity F{u) = uu^ is locally Lipschitz continu- 
ous from L'^{K) H~"' for some < 7 < 1 . Then by a cutoff technique we 
have the existence and uniqueness of local solution for some stopping time T 
by the classical application of Banach Fixed Point Theorem for stochastic 
evolutionary equation [37]. The global existence and uniqueness is thus fol- 
lowed by the a priori estimates in following part. □ 

To construct a stationary solution by Bogolyubov-Krylov method, we 
need the tightness of {C{u{T,^))}r>o in the space L'^{K). By the com- 
pact embedding of H^{K) C L^{K), we just give the uniform estimates 
of {m(t, ^)}t->o in the space H^{K). For this we need some estimate in the 
spaces L°°(M) and L^iK). 

2.2.1 Estimates in the space L°°(]R) 

We adapt the approach used by Zuazua [HUj for a scalar convection-diffusion 



equations to give an L°° (M) estimate for solution to equation ( 10 ) with initial 
value u° G L'^{K) f] L^{R). 
For this we introduce 

^ fl, M > 0, Jl, M < 0, 

sgn(M)+ = <^ and sgn(M) = <^ 

I 0, M < ; I 0, M > . 

Then for u G ^^(IR) with Vu(t) G iv^(]R), J^U!:^(p{u) d^ < for any non- 
decreasing {p G C"'^(]R). By a density discussion J^u^^sga^u)"^ d^ < 0. 
Moreover, J^uu^sgn^u)"^ dC^ = and J^{C,u^ + u)sgn{u)'^ dC^ = 0. Denote 
by = sgn{u)^u. Now let m(r) = ||m°||l°°{r) + lh(''")||L°°(R) with ?7(r) 



solves (14) with a = . Then multiplying sgn(-u — m)+ on both sides of (10) 



and integrating on M x [0, r] with r > 

{u{T,O-m)+d^<0. 
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Therefore u{t,C,) < m(r) for any r > 0. Similarly u{t,C,) > —m{T) , r > 0. 
Then ||u(r) Hiooj^K) < m(r) , for all r > 0. Moreover, by the property of r/°(r) , 
Em(r) is bounded by some positive constant. 



2.2.2 Estimates in the space H^{K) 

We first give estimates of solution {m(t, ^)},->o in the space L'^{K). 

First for u{t,^) solving (10) with m° G L°°(]R) D L'^{K) we decompose 



u{t, ^) = uo{t, ^) + Us{t,^) with Mo G Eq and Ug & Eg ■ Then by the assump- 



tion (11) 



duQ 



[Cug- Ps{uu^)]dt + dW: 



2 ■ 



Then 



uo{t,0 = uo{0,0 = (n(0),eo)eo(0 



So uo{t,C,) is totally determined by the mass of initial value which is denoted 

Now introduce v = u — 7]°' , then we have the following random evolution- 
ary equation 

Vr = Cv-{v + 7]''){v + 7]'')^ + ar]'^. (16) 

We first give a uniform estimate on v{t) in the space L'^{K). Similarly write 
V = Vq + Vg E Eq (B Eg , then 

vo{T,0=Mr,0 = Meo{0- 
Multiplying v in the space L'^{K) on both sides of the equation (16), 

Consider the third term on left hand side of the above equality 



d 



{v + ry-{vK)d^ 



- I v'^v^Kdi- I vri'^v^Kd^ 
Jr Jr 

- [ v^ri^Kd^- [ rj'^ri^vKd^. 
Jr Jr 



We estimate the left four terms separately. By Cauchy inequality for any 
e > , there is positive constant such that 



vrj^v^K dC, 



< \\v'^\\l^\Hlhk)\\v^\\l'Hk) 
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Jr 



< \\v\\L-°\\v\\L2(K)\\'n'^\\mK) 

< \\v'"\\miK)\\v\\l2^K), 



\l^{k) 



< 



< a 



l|4 



^\\l2(K)\\v\\l2(K) 
2 



+ e V 



Integrating by parts 

[ v\Kd^ 
Jr 



v^Kdv = -2 I v^v^Kd^- [ v^K^d^. 
Jm. Jr 

By property |6] in Lemma [T| for any e, e' > 0, g > 2, there exist positive 



constants , C^'^q and R such that 



[ v^K^di = - [ viK^'^v^K^/^di 
Jr 2 J^ 



< 



v^eKdi 



- 1/2 






/ v^Kdi 
Jr 



1/2 



< C\\v4l2(k)\\v\\l^{k)\\v\\l--{r) 

< ^^CMhiK) + [^'MhiK) + a',,l|t;||L(i,(o,«))] \\v\\ 

< 3 [eC + e'CJf |||,x=(]R)] ||f|||i2(;^) + 3CeCe',g,ij||f lli^(iR) 

with some positive constant C^i^q^R , then 



f v\Kd^ 
Jr 



\ [ v'K^di 

•-> JR 

< [tC + e'a||t;||2^(jjj ||t;^||i2(x) + aC,',g,fl||w||io 



Then for any e and e' > , there is a positive constant that we still denote 
by Ce such that 

+ ^J^"lli-(R)ll^lli2(if) + h'"\\H\K)\Hh(^K) + e||^^lli2(^) + 
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+ \\v''\\hhk)] WvsWl^^K) + [-'2 + <1 + C) + e'ail^llioo(K)] \\vo\\m(K) 



+ [26 + ah"llio.(M) + hiiHi(x)] 

\m{K) + "C'e||?7"||^2(;^) 



L2(X) 



Now choosing e and e' > small enough and noticing that ||f < 
2|klli?i(E-) , 

where 

ff(r,a;) = 2 [|(3 + C) + e'a||t;||i.(i,) + + Hr/H^Hx)] 

and 

Mr,u;) = [-l + H{T,u)]M' + aC,\\r{T)\\l,^^^ 

\h^{K) 

Then by the Gronwall inequality for any r > 



+ C.\\v''{r)\\%,iK)+CeC,,,MHU 



\vs{r)\\l,^K) < e'2^+/<:^W'^^||^;(0)||i.(^) (17) 



+ 



By the construction of ri°'{T,uj) and the estimate on u{t) in L°°(]R) (section 



2.2.1) , choose e and e' > small enough there is a random variable tq{u) 



such that for r > tq almost surely 

exp |-|r + ^ H{s) rfsj < exp {-|r} . 

Then v{t) is time uniformly bounded by a tempered random variable Ri{6t-u)) 
in the space L'^{K). 

Remark 9. Notice that 1]°' is very small in L^{Qq, H^{K)) by choosing large 
a > (Lemma |5|. Then by the definition of h and estimates of u{t,^) in 
the space L°*^(M), for small initial value uq G L'^{K) fl L°°(]R), h is small 
in L^(f2o)- Then by estimate (17) for r > tq is large, v{t) is small in the 
space L'^{K) almost surely. 

Now given any ri > , in the mild sense 

v{t + ti) = S{t)v{ti) + ^ / S{t - a)v{Ti + a) da 

Jo 
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Jo 

X {v{ti + cr) + t]'^{ti + a))/: da 

+ a S{t — a)r]°'{a) da . 
Jo 

Then taking (K) norm in the above equation, by the growth of the semi- 
group, for some positive constant C 

\\v{r + Ti)\\Hi(K) 

< C{l + ^)\\v{n)\\mK)+ r {l + ^^)e~^'^Ri{e^uj)da 
A Jo ^T-a 



+ C [ {l + —^)e'^/^Ri{e^uj)\\v{n + a)\\H^K)da 
Jo Vr - (T 

+ C [ (1+ / )e-^/''{l + Riie^uj))r]''{e^u)da. 
Jo \Jt -a 

By Gronwall lemma, and the tempered property of Ri{6rU}), there is a tem- 
pered random variable R2{0rU}) such that 



for any Ti > . Then we have the uniform estimate of v{t) in the space (K) 
and the compact embedding of H^{K) C L'^{K) with the property of 7]°" 
yields the tightness of {C{u{T))}r>o^ the laws of u(t), in the space L'^{K). 
Then the classical Bogolyubov-Krylov method [1] yields the existence of an 
invariant measure denoted by /i. 

Now choose random variable with C{u^) = fi , then the solution u{^, u) 



with initial value m° is a stationary solution to (10). By the self similar 
transformation for t > 1 , 

u(t,x,uj) = -j=.u (^x/Vi,u?j = — (^x/Vi,u?j 

is a random self-similar solution to stochastic Burgers' equation ([5]). 
Remark 10. By the stationary property of u, ^/iu{t,x,uJ) just depends on 
^ = x/y/i in the sense of distribution. So u{t,x,uj) is called the random 
self-similar solution to stochastic Burgers' equation 

2.3 Locally asymptotic convergence to self-similar so- 
lutions 

Self-similar solutions are important in study the dynamics of the system (|5]). 



Next we show that solution u{t, ^) to equation ( 10 ) tends to a unique station 
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ary solution u{t, as r — )■ oo under some conditions. This shows the asymp- 
totical convergence of stochastic Burgers' equation ^ to the self-similar so- 
lution. 

First we have the following result on the stationary solution to equa- 



tion (10). 



Lemma 11. Any stationary solution, which is small in L'^{Qq, L°°{M.) fl 
L'^{K)), to (10) is uniquely determined by it's mass, that is for any given 



M G M , there is a unique stationary solution u{t,^) to (10) with 



Proof. Suppose Ui and U2 are two stationary solutions to (10) with 



Jm. Jm. 



(19) 



Let U = Ui — U2 , then 

Ur = CU- (nif/)^ - UU^ 
By the Cole-Hopf transformation 



we have 



y = f/exp|-i f/(y)rfy| 



(20) 



(21) 



Notice (19), V has zero projection to Eq, that is V ^ Eg . Multiplying V on 



both sides of (21) in L\K) 



^ d 



w 



1||T/||2 



^11 " Wl^k) - -2\\ " \\m(K) 



[ uiVV^Kd^-^ I uiV^^Kd^. 
Jr Jm. 



For the last two terms by Cauchy inequality for any e > , for some positive 
constant 



uiVV^Kdi 



< CA\ui\ 



L°°( 



and 



/ 

JR 



uiV'^iKdi 



< II^iIIl-(m) / V^iKd^ 
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< ||^^i||l°o(r) 

< C'ell'l^l ||ioo(• 



\^\\h{K) + ^Il^clli2(ii-) • 



Then by ||V||l2(;^) < 2\\V\\hi(k), 
d 



for small e 

(-| + |e + 3CJni||i. 



Now for small Ui in L'^{VLq, L°°(]R) fl L'^{K)), by the stationary property and 
the Gronwall inequality, almost surely 

V{t) , r oo 



which yields the uniqueness of the stationary solution satisfying (18). The 
proof is complete. □ 

By the above result we show that the long time behavior of some solution 
of ( 10 ) is approximated by a unique stationary solution. For this we first show 
that the stationary solution u constructed by Bogolyubov-Krylov method is 
bounded by the bound of u{t,C,). 

Lemma 12. For any solution u{t,C,) to equation (10) with initial value vP G 
L^{K) n L°°(M), and K\\u{T,^)\\'j^2(K)f^x.°°(^^-j < C for some t -independent 
positive constant C , then there is a stationary solution u such that 

u{^)d^= u{T,^)d^ and ^\\u\\l2^K)nL--{R) < C . 
Jm. Jr 

Proof. Denote by /i^ = Cu{s) the distribution of u{s, ■) in the space L'^{K). 
The Bogolyubov-Krylov method introduces the following probability mea- 
sure on L^{K), 

1 



1^7 



jjs ds . 



and finds a limit point of {/2^}^>o • By the estimates of solution in the 
space H^{K), {/2^}^>o is tight in the space L'^{K), then there is a probability 
measure fi on the space L'^{K) and subsequence r„ with r„ — oo , n — t- oo , 
such that 



in weak sense [5] . Now let u be the stationary solution to equation ( 10 ) with 
initial distribution fl, then for for any r > , 



E ||M||i2(x) 



\u{^,^;'v)\\h{K)Kdv) 
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wherever u{-,^]v) is the solution to equation (10) with initial value v G 
L'^[K). By the construction of /i 

^\M\^(K) = lim — /" f \\u{T,i]v)\\l2,K^^^s{dv)ds 

1 



n— >oo 7", 



Tn 



- ~ / IE||u(r,^;M(s,0)|li2(x)C?s 



n JO 



< c. 



Similar E ||m||^oo(]r) < C . The proof is complete. □ 

Now we prove the following local asymptotical convergence to self-similar 
solution. 

Theorem 13. For any solution u{t,C,), which is small in L'^{Qq, L°°(R) fl 
L'^{K)) for any t > , there is a unique stationary solution u{t,^) such that 
almost surely 

||n(r) - n(r)||j:2(^) ^ 0, r ^ oo . 

Proof. Assume for any r > 0, E ||M(r, ,^)||^2(;^) < C for some positive con- 
stant C. Let 

M = [ n(r,Orfe, 



then by Lemma [TT] and Lemma 12 , there is a unique stationary solution, 
denoted by u{t,^), satisfies 



Moreover, u is small in the space L'^{^Iq, L'^{K) fl L°°(]R)). 
Let U (r, = u{t, ^ - 0> then 

Ur = CU - {uU)i: - UU^ . 



By the Cole-Hopf transformation (20), we similarly have 

Vr = CV- (uV)^ . (22) 
By the choice ofUjVEEg, and by same discussion in the proof of Lemma 11 , 

^II^IIl2{A-) = (~i + + 3Ce||u||^oo(]f5)) ||V^||i2(;^) . 

By the stationary property of u{t,C,) almost surely 

||F(r,OIU2(K)^0, r^oo. 

The proof is complete. 

□ 
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Then for stochastic Burgers' equation (|5| on time interval t > 1 
Theorem 14. For t > 1 , if the solution u(t, x) to equation ^ satisfies that 
e^/2u(e^,e^/2^) zs small zn the space L'^(yto,L\K) n L°°(M)) for any t > , 
then there is a unique stochastic self-similar solution u such that for almost 
all u E fl 

Vi\\u{t,X,Uj) - u{t,X,Uj)\\L2(R) , t ^ oo 

where ^ 

u{t,x,uj) = —=u (x/Vi,io 



with u is the unique stationary solution to equation (10) with J^u{^) = 

Remark 15. By the above discussion, if u is smaller in L'^{Qo, L°° {R)nL'^ (K)) , 
e can be chosen smaller and then the exponential converge rate is closer 
to -1/2. 

2.4 Asymptotic convergence described by a local ran- 
dom invariant manifold 

Next we show the locally asymptotical convergence of the solution u{t, ^) by 
the random invariant manifold theory [H UHl e.g.]. This approach does not 
require the solution u{t,^) is small in the space L'^{K) fl L°°(]R). 

Notice that the nonlinearity uu^ is local Lipschtiz, a global random invari- 
ant manifold is difficult to be constructed [iHl HE] . Recent work by Blomker 
and Wang [8] gave a cut-off method to construct a local random invariant 
manifold (lrim) for SPDEs with quadratic nonlinearity. Here we also consider 



a LRIM for equation (10) by a cut-off technique. 



We construct a LRIM for a stationary solution u . Notice the noise is 



additive in (10), for any stationary solution u we introduce 

Then denote by B{u,v) = ^{uv^ + vu^), 

Vr = Cv - [B{v + u,v + u) - B{u,u)] (23) 
which is a random evolutionary equation. By the property of C, deterministic 



approach for the well-posedness yields that equation (23 ) defines a continuous 
random dynamical system (p{T,u) with driven system {6*^},-. 

Now let Pg be the projection from L'^{K) to the space Eg, and 

V = Vq + Vs e Eq ® Es . 
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Then 



t'o = 

Vs = Cvs — [B{v + u,v + u) — B{u,u)] . 

Notice that nonhnearity 

F{v, OrUj) = B{v + u,v + u{t)) - B{u, u) 

is non-Lipschtiz, to construct a LRIM for ip{t,uj) we need some cutoff tech- 
nique [8j . Denote by i?/j(0) the ball with radius R that is 

Br{Q) = {u e L\K) : \\v\\l2^k) < R} . 
For any set A C L'^{K), define the following distance 

dist(M, A) = inf ||n — f 

for any u G LP'{K). Then introduce the following cutoff function on the space 
L\K) 

B^{u, u) = Xr{u)B{u, n), ue L\K) 

where i? > and Xni!^) = x('"/-R) with x '■ L^{K) — > M is a smooth bounded 
function with xi'U') = 1 if ||'u||i,2(^-) < 1 and x('u) = if ||m||2,2(^) > 2. Then 
the nonlinear term B^{u,u) : L'^{K) — )■ H^'^{K), < 7 < 1 , is Lipschitz 
continuous. More precisely for some positive constant Cb 

\\B{ui,Ui) - B{u2,U2)\\h-<{K) < 2RCb\\ui - U2\\l2{K) 

for all ui,U2 G L'^{K) with ||^ii||L2(x) ^ R and ||^^2||L2(_ft:) < R- Now define 
F^{v, OrUj) = B^{v + u,v + u)- B^{u, u), 

then 

^^PLHK),H-'{K)iF^iu,u)) = Lr := 2RCb ■ 
Consider now the following cut-off system 

= (24) 
vf = Cvf - F^{v,e^u) (25) 

which defines a continuous random dynamical system <f^{T, u). By the spec- 
trum property of C, for small R > , the Lyapunov-Perron method for SPDEs 
is applicable to the cutoff system (24)-(25) [HI [iHl e.g.] . Then for enough 
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small R, (p^{t,uj) has a random invariant manifold A^^tl'^) which can be 
represented by 

■^cuti^) = {(^0, h{vo, u)) : vo e Eo} . 
Here h{-,u) : Eq ^ Eg is Lipschitz continuous with h{0,u) = and 



h{vo, co) = - e-^'F'^iv; + vl, e,u) ds 

J ~oo 



where {vq^v*) is the unique solution of 
fo(r) = Vq 

J —oo 

in the Banach space 



(26) 
(27) 



C; 



V e C((-oo,0],L2(K)) : supe"^^||t;(r)||i2(j^) < oo 

T<0 



endowed with norm 



sup e 

r<0 



-An 



v[t)\\l\k) , 



where — 1/2<A<0. In fact for any vq G Eq define nonlinear operator 
r:C,-^C,-by 



(wo,0 ^ riv*Q,vl;vo) = right-hand side of (26)-(27) 



Then for i? > is small enough, a direct calculation yields the mapping T 
is contraction; that is, T has a unique fixed point {vq,v*) = {vo,v*) G 



which is the unique solution to (26)-(27). 
Then 



MU^) = A^Ko;) n 5^(0) 



(28) 

defines a LRIM for ip{t, u). Furthermore, by the similar discussion for stochas- 
tic Burgers' equation on bounded domain [8] , the random invariant mani- 
fold A^^^.(ci;) is almost surely complete. That is, we have the following result. 
Theorem 16. Assume R > is small enough. Then for any solution 
'^'^('^jO = ("^(f ('^7 0' "^f (''"i 0) of the cutoff system (24}-{25) there is one 
orbit V^{t,^) on Ai^ut ^^^^ K)^ = "^o^ ^^(^^ ^^^^ 

IIKV),^f (r)) - (nV),KV))llL^w 

< ||t;^(0) -\/^(0)||e-^*" 
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where 



X* 



4' 



1 + 



for some S > . 

The above result yields that 
Theorem 17. For any vP G -B_r(0), 



dist(v9(r, A^ioc(^T^)) < 2i?e 
/or a// r < ro(a;) = inf{r > : ip{T,u)uQ ^ Br{^)} . 



-A*7 



Theorem [16] and [17] describe the local attractive property of the zero so- 
lution of (23). Especially for the solution (0,fs(r, ,^)) to (23), by h{^,uj) = , 
the attractive orbit on invariant manifold A^ioc(i^) is the zero solution. That 



is the solution u{t,^) = {u ,Us{t,^)) to (10) is attracted by a stationary 
solution u{9rUJ,0 = {u^,Us{9rUJ,C)) for r < ro(ci;) provide m(0,^) — ^(0,^) 
lies in a small ball -B_r(0) C L'^{K). By the construction of the local random 
invariant manifold, u is the unique stationary solution that attracts u. We 



then still have the result of Theorem 13 in a small ball of L'^{K) fl L°°(]R). 
However, the local random invariant manifold method does not restrict so- 
lution u to be small in the space L'^{K) fl L°°(]R). We draw this conclusion 
in the following corollary. 
Corollary 18. For any solution u{t,C,) 

let u{t,^) 



u 



[t,^)) to equation (10) with 
{uQ,Us{9rU!,^)) be one sta- 



tionary solution of equation [T^ - Then ifu{Q, ^) — u{0,^) lies in a small ball 
-B_r(0) C L'^{K) with R small enough, almost surely 

\\uiT,^)-uiT,mLHK) < \H0,O-u{0,mLHK)e-^''^ 

for r < To . 

Now the following theorem applies to the stochastic Burgers' equation dS). 
Theorem 19. For any solution u{t,x) to stochastic Burgers' equation (pj, 
let u{t,x) be a self-similar solution to equation Then ifu{l,x) — u(l,x) 
lies in a small ball -B/j(0) C L'^{K) with R small enough, almost surely 



\/i\\u{t,x) — u{t,x) 







t — i- oo . 



Remark 20. The local random invariant manifold we constructed here is 
a global one for the random dynamical system ip{t,u) by the next global 
asymptotic convergence. 
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2.5 Globally asymptotic convergence to self-similar so- 
lutions 



We show the asymptotic convergence in probabihty of any solution of stochas- 
tic Burgers' equation. For this we first consider the Markov semigroup defined 



by the solution of the equation (10). 

Denote by Ai the space consisting all probability measures on space 
L"^ [K) n L°° [M.) and endow Ai with the topology of weak convergence. Define 
continuous Markov semigroup {Pt}t>o on Ai 

Pri^iA) = P{u(r, ■)eA}, fiEM 

for any Borel measurable set A C L'^{K) fl L°°(]R) and u(t, ■) is the solution 



to equation (10) with initial value vP distributes as // . The space Ai is too 
larger for our purpose. For this we introduce the following subspace 

Ai2 = l^i E Ai : I \\'^\\'L■^(K)^^{du) < oo > . 



By the discussion in section [272] , there is a /i G A^2 such that Prp- = fi , 
which is called stationary measure of Pr . We next show that for any yU G A^2 
there is a unique stationary measure fl G Ai2 such that PrfJ. converges weakly 
to /i as r — 7- oo , that is 



f{u)Pr^{du) — )■ / f{u)iJ,{du), r — !■ oo 

L^{K)r\L°°{R) J L^{K)nL°°{R) 

for any bounded continuous function / : L'^{K) fl L°°(]R) — )■ M . 

Associate the solution to equation ([To| we choose fj. G Ai 2 which has the 
form 

fi = 6m* fJ-s (29) 

where Sm is some Dirac measure on Eq and fig is supported on Eg . Then 
consider the limit of Pr^i as r — t- 00 . First by the estimates in section 2.2 we 



have a measure /i and subsequence r„ with r„ — )■ 00 , n — )■ 00, such that 

P^-nA^ 1 n — )■ 00. (30) 



To show the uniqueness of /z need a contraction property of the system (10). 
This is from the deterministic result . 
Lemma 21. For any , v? e L'^{K) fl L~(M) w%th 



Jr 
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Let u^{t,^) and u^{t,^) be the solutions to equation (10) with initial value 
and respectively. Then the function 



/ \u\r,i)-u\r,i)\di 



is strictly decreasing almost surely. 

Proof. Let U{t,C,) = u^{t,^) — M^(r, then we have the following linear 
equation 



u 



mi- 



By the estimates in section 



2.2 



[u^ + V? — is bounded by a random 
constant. Then for any a fixed oj &VLq , the result is followed by the discussion 
for deterministic system [50l [30] . □ 



Now for any ^ & M. with form (29 ), let /i be a stationary measure and (30 ) 



holds. Suppose /i' is another stationary measure of P.^ such that for some 



Pr' fi ^ fi' , n — )■ oo . 



(31) 



Denote by u(r, ^) and u'{t,^) the solutions of equation ( 10 ) with initial value 



u^{C,) and u'^iC,) distributes as fi and fi' respectively. Then 



By Lemma 21 , the function 



[ |w(r,0-^'(r,OMe 
Jr 



is almost surely strictly decreasing in r which contradicts stationary of u 
and u'. Then we deduce the following result. 

Theorem 22. For any vP E L'^{K) fl L°°(]R), the solution u{t,^) to equa- 
tion (10) with initial value converges in distribution , as r oo , to u in 



space L {K) which is the unique stationary solution to equation (10) with 



Jr Jr 



Next we show the above convergence in distribution is in fact a conver- 
gence in probability. For this we need the following fact stated by Gyongy and 
Krylov [23]. Let X be a Polish space with Borel sigma algebra. A sequence 
{Xn} of X- valued random variables converges in probability if and only if for 
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every pair of subsequences {X^} and {Xi} , there exists X x X- valued sub- 
sequence Zk = (Xm(fc), converging in distribution to a random variable 
Z supported on {{x,y) E X x X : x = y} . 



Now for solution u{t,^) to equation (10), consider any subsequences 



{ui'^ni)} and {u{ti)}. By the Theorem 22, both subsequences converges in 
distribution to a random variable supported on 

{iu,v) e (L^K) n L^(R)) X {L^{K) n L°"(M)) ■.u = v}. 

Then we draw the following result 

Corollary 23. For any solution u{t,^) to equation (10) with initial value 
vP G L^{K) nL°°(]R), there is a unique stationary solution u{t,C,) such that 

\\u{t,^) — u{t,^)\\i2(^j^-j — i- , in probability as t ^ oo . 



Remark 24. By the random invariant manifold discussion in section 2.4 
the convergence rate approximates 1/2 after a long time. Furthermore, by 
the above global asymptotic convergence result, any stationary solution is 
uniquely determined by its part in Eq, that is for any uq G Eq , there is 
Us = h{uQ,uj) such that {uQ,h{uQ,uj)) is the unique stationary solution with 
Uq part in Eq . Then by the local random invariant manifold discussion in 



section 2.4, equation (10) has a global random invariant manifold which can 



be represented by {(mq, h{uo,oo)) : uq G Eq} . 

Then the following theorem applies to the stochastic Burgers' equation pi) 
Theorem 25. For any solution u{t,x) to stochastic Burgers' equation 
there is a unique self-similar solution u(t, x) to equation such that 



Vi\\u{t,x) 



u(t, x) 



, t ^ oo , 



in probability. 



3 Self-similarity emerges in examples of prac- 
tical interest 

The proof of the previous section [2] shows one case where we can prove the 
emergence of stochastic self-similarity. However, the stochastic slow mani- 



fold framework of section 1A_ strongly indicates that stochastic self-similarity 
emerges in a much wider class of stochastic systems. We proceed in the re- 
maining sections to explore, albeit less rigorously, several example spdes of 
the form (|6|) and of significant physical interest. 
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We investigate particular stochastic systems and the self-similarity in 
their stochastic slow manifolds by writing the dynamics in the Hermite ba- 
sis: from the eigenbasis ^ we express the transformed field as the spectral 
expansion 

oo oo 

u{t,0 = J2uk{T)e,{0 = Y,u,{T)ckH,{^/V2)G{0, (32) 

fc=0 k=0 

for normalisation coefficients Ck and mode amplitudes Mfc('r). We also write 
the Wiener noise in the cylindrical expansion Q. The very first differ- 
ence with the rigorous analysis of section [2] is that in this section we allow 
the noise to have a mean component in space; that is, in this section the 
noise coefficient bo . The consequent direct forcing of the fundamental 
Gaussian mode makes an immediate qualitative difference to the long term 
evolution that applied scientists and engineers will appreciate. 



3.1 Stochastically forced diffusion 

Obtain the simplest stochastic self similarity for additive noise and when 
there is no nonlinearity: in this subsection we set / = and g = 1/t in 
SPDE ([!]). Substitute the Hermite expansions (32) and (|9| into the trans- 



formed SPDE (|6]) and equate coefficients of the basis functions ek{^) to find 
the decoupled SDE system 

ilk = -^Uk + bkWk, /c = 0, 1,2, . . . . (33) 



For mode numbers k > 1, the SDE (33) describes an Ornstein-Uhlenbeck 
process and so the amplitudes of the fast modes are Uk = Mfc(0)e~^'^/^ + 
hkG~^^^'^-^Wk where we define the stochastic convolution 



which is the Ornstein-Uhlenbeck process satisfying the SDE dz = —(5z dr + 
dwkir). Consequently, exponentially quickly Uk — )• hkC'^'^^'^-kvok as log-time 
r — )■ oo , and this exponential approach is 0{e~'^^'^^ due to the decay of the 
leading 'fast' component of the initial compact release. 

But the slow mode satisfies Mq = ^qUIq with solution uq = aij) = boWo^r). 
Hence, in this case the emergent stochastic slow manifold of the spde (j6]) is 

oo 

u = -^e-^'/' + y2ek{0e-'-^'^w. 
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where the amphtude a = boWoir) = boWo(\ogt). That is, under the direct 
forcing of the Gaussian structure, the stretched field u undergoes a random 
walk in the amplitude of the Gaussian, while exhibiting zero-mean fluctua- 
tions due to the past history of the other noise components. Since the Wiener 
process wo{t) = (9(r^/^), this predicts the original stochastic self similarity 
field will be u = C(ri/2(logt)i/2) as t ^ oo . 

For comparison, recall that section |2] proves the emergence of a stochas- 
tically stationary distribution in the case when there is no mean noise com- 
ponent. 

The self-similar random walk emerges from transients of relative magni- 
tude C(e-^/2^ = 0{t-^/^); that is, of absolute magnitude in field u of C(t-i). 
Section |4] returns to this case to argue that stochastically moving the space 
origin and stretching time empowers us to eliminate the leading two stable 
modes ui and U2- The new view of Section|4]is a stochastic self-similarity that 
emerges somewhat quicker, with relative transients (9(e~'^'^/^) = 0(t~^^'^y, 
that is, of absolute magnitude in field u of 



3.2 Cubic reaction enhances decay 

Here consider the case of cubic reaction, / = — u^, and additive noise, g = 
1/t . The stochastic slow manifold shows that not only does the cubic reaction 
aid the decay, but also noise-noise interactions increase the exponent in the 
self-similar decay. 



3.2.1 Change to Hermite basis 

For simplicity, initially just project the dynamics of the spde ^ onto the first 



three Hermite modes; section [3.2.3| uses computer algebra to implement more 
modes. Recall the Hermite functions are -fffe(C) = (— l)'^e'» /'^d^e~'^ 1'^ jdC,^: 
the first few are -f^o(C) = 1 ; -f^i(C) = C and -^2(0 = — 1 • As given by 



equation (32), we expand the solution field in the corresponding basis. The 



complication is the cubic reaction term 

£-3^2/4 

l,m,n 



We want to expand = J2k dk^kiOy take the weighted inner product of 
this cubic reaction with = CkHk{^)G{^) to determine the coefficient 

dk=22 '^lUmUnCkClCmCn J H k{^) H H ^{^) H ^^^^ 
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l,m,n °° 

upon substituting ^ = a/2/3 C • Expressing the Hermite functions as polyno- 
mials in (, expanding the quadruple products into a high order polynomial 
in C, and then recasting the polynomial in a sum of Hermite functions, the co- 
efficient of the Ho{() mode determines the above integral. Omitting details, 
upon truncating the above sums over modes to just the first three modes, 
k, /, m, n = 0, 1, 2 , we find coefficients in the cubic are 

do = -^{^ul + luoul - ^ul + \uou\ - ^^0^2), 

cii = -^{\ulu^ + \u\ + 

VOTT 

d2 = ( - + \ulu2 - "fuoul + l^ul + \ulu2) . 

V ovr 



Consequently, the stochastic system, when projected onto the first three 
modes, is approximated by the set of SDEs 

Uq= + boWo + ^= ( - ^Mq - ^UquI + - \uquI + -^ulu2) , 



ill = -^Ui + biWi + -7={- ^ului - - luiuj), 

V OTT 

U2 = -U2 + + ^=(^^0 - i«0"2 + ^UquI - - |miM2) . (34) 

Computer algebra checks this derivation, and also computes the correspond- 
ing systems for higher order projections. 

3.2.2 A normal form separates fast and slow modes 

As a preliminary to the more complete construction of the stochastic slow 



manifold in section 3.2.3, this subsection shows that a stochastic coordinate 



transform [H Ch. 8] separates the stochastic slow and fast modes in the 



projected system (34). 



In the projected system (34) the dynamical variables ^^(t) are linearly 
diagonalised with constant coefficients, and the nonlinearities are of multino- 
mial form. In such a case, constructing the necessary stochastic coordinate 



transform is routine [H]. Indeed, a web service [12] analyses the system (34) 
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to construct a near identity, stochastic, coordinate transformation from vari- 
ables Ukir) to new variables Ukir). The constructed transformation, begin- 
ning 



Ui 
U2 



Von 



(35) 



transforms the system of SDEs (34) to the equivalent system 



3^Ui - \UlUi 



V2bob 



2Woe 'kW2 



371U2 



- - lU^U2 + ( - lU^b2W2 + ^UoUib2W2) 

+ U2 [-\b\wie^''^'^'kwi - Y^blw2e~'^icW2 + 1 60^2^06"'^ ^1^2] • (36) 



Alternatively, one can straightforwardly confirm the order of accuracy of 



SDES (34) 



(35) and (36) by simply substituting them into the projected system of 



The long term dynamics are readily apparent from the transformed sys- 
tem (36). We immediately deduce the existence of a slow manifold, its emer- 
gence, and its evolution. Observe that Ui 



is invariant in the 



transformed system (36); a stochastic coordinate transform such as (\35h may 



always be found to create such invariance to any specified order [21 [271 E] . 



Due to the exponential decay of Ui and U2 in the deterministic parts of (36), 
provided the magnitudes b^ of the stochastic effects are not too large, the 
invariant manifold Ui = U2 = will be almost surely exponentially quickly 
attractive. The emergent stochastic slow manifold is thus Ui = U2 = . 



Evolution on the stochastic slow manifold is given by the first line of (36), 
for all time. Substituting Ui = U2 = in the transform (35) then gives the 
shape of the stochastic slow manifold in the w^-variables. 



The transformed SDEs (36) indicate that the cubic nonlinearity enhances 



the rate of attraction to the stochastic slow manifold: the deterministic part 
of the two fast modes are Uk ~ + ^Uq / {\/3Ti)~\Uk ■ But the most 

important aspect is that, by continuity, there exists a finite domain near the 
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u = equilibrium such that only very rare stochastic events could overcome 
the exponential attraction to the stochastic slow manifold. Thus we expect 
the stochastic slow manifold to almost always emerge from some finite domain 
of initial conditions. 

Consequently, the following 'shadow' modelling applies. The SDE sys- 



tem (36) and the stochastic transform (35) together describe a stochastic 



process in the Ukir) coefficients of the Hermite basis functions efc(^), and 
thus describes a stochastic process in the field u{t,^). By the asymptotic 
construction, this stochastic process is in an asymptotic sense 'close to' or 
'shadows' the original spde (|6]), especially when carried out to better reso- 
lution as we do in the next section |3.2.3[ Now, in physical applications the 



original coefficients, functional forms and noise spectrum in the spde ^ 
are never known exactly. Thus in physical applications, predictions deduced 



from systematic shadowing stochastic processes such as (|36|) and (|35|), are as 
useful as predictions from the spde 



On the emergent stochastic slow manifold, the system (34) evolves in the 



long term according to the first SDE of the system (36). The most important 



part of the quadratic noise-noise interaction terms in (|36|) is their effect upon 
the mean drift. 



Analysis of such noise-noise interactions [ID] shows they 
generate drift and fiuctuations: neglecting these fiuctuations the slow mode 



of the SDEs (36) becomes 



(37) 



where we use Uq = a to denote the amplitude of the stochastic self similar 
solution. In this SDE: the direct forcing hoWQ promotes a random walk among 
the self-similar profiles, as for linear diffusion; the cubic reaction — \a^ reflects 
the cubic reaction of the original physics; but the —{\b\ + \h'^a term accounts 
for noise-noise fluctuations enhancing the exponent of the similarity decay 



rate. Just solving -\/37rd ~ + l^i)*^ gives amplitude a oc e "'^ = t " 

= {^bf + \b2)/\/S7r and hence predicts the original fleld 
-1/2-0 Physically, because of the nonlinear shape of the 



for exponent a 
u oc t~^/'^a oc t~ 

cubic reaction, and in comparing fluctuations that enhance the local fleld 
with those that decrease the local fleld, the flrst generates reactions that are 
slightly larger than the reaction is decreased by the second. Thus such a cubic 
reaction enhances the similarity decay rate through noise-noise interactions. 
Our systematic resolution of the noise-noise interactions discerns this noise 
enhanced decay. 
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3.2.3 The emergent slow manifold of stochastic self-similarity 

When, as for the spde ([6]), there is a large number of noise excited fast modes, 
then the full normal form coordinate transform is impossible to construct. 
Fortunately, we just need to construct the stochastic slow manifold part of the 
coordinate transform, as is summarised in this section. Here we project the 
SPDE ^ onto the first nine modes and construct and interpret the resulting 
stochastic slow model. 

Computer algebra obtains a system of SDEs for Uq, . . . ,Us, analogous to 



the system (34), but which are far too involved to record here. Established 
methods |10l E] , available in computer algebra via the web [33] , then anal- 
yses the system of SDEs to determine that the stochastic slow manifold is 
approximately 



/37T 



1 — t/2 • ^ 1 — T 

Ml ~ bit ' -kWi , U2 ~ H ^= + 026 'kW2 , 

3v67r 

^3 

M3 ~ b-^e~^'^^'^'kW3 , U4 ~ ^= + b4^e~'^'^'kW4 

18v 27r 



3 



Vs. 



3 



M5 ^ 656 ^''^'^-kw^ , Uq ^ H -1= + b^e ^^-kWe - 

81v37r 



7 -7t/2 . y/TOcL^ _A 

o48vo7r 



As shown explicitly for the low order system in section 3.2.2 by continuity, 
and except for rare events, this stochastic slow manifold will almost always 
emerge from all initial conditions in its neighbourhood. 

The methodology implemented in the web service [13] also constructs 
the stochastic evolution on this stochastic slow manifold. Effects quadratic 
in noise are very complicated — too complicated to record here — due to the 
need to resolve the multitude of noise-noise interactions that occur in the 
fast modes [IQ]. Here we retain only their cumulative drift effects of the 
noise-noise interactions. Computer algebra then finds the evolution to be 



the following more complete version of the earlier (37): 
\/37rd ~ — \a? + \/37rboWQ 
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a LgOl -h + 54^3 + 216^4 + 270"5 + 972*^6 + 3402"7 + 1296"8j • 

(38) 



The first line of tlie stocliastic slow mode (38) contains the direct effects of 
the cubic reaction and the stochastic forcing. The second line is a multi- 
plicative noise term that could be replaced by one independent noise term 
with volatility coefficient being + 61/24 + 56^/729 + 7061/46656. The 



last line of (38) enhances the self-similarity decay rate through noise-noise 
interactions. 

A Domb-Sykes plot [3^ e.g.] of the ratio of the coefficients in the last 
line suggests the corresponding infinite series converges provided assump- 



tion (13) holds. However, extant stochastic slow manifold theory is limited 
to effectively finite dimensional dynamics such as the spde ([6]) projected onto 
the first nine Hermite modes that we analyse here. Nonetheless, in princi- 
ple we could construct the stochastic slow manifold model to some level of 
approximation for any finite truncation of the noise ([9]). 

3.3 Modelling other stochastic systems 



Section LI established that a stochastic slow manifold approach illuminates 
a broad class of stochastic reaction-diffusion pdes. This approach is then 
supported by some rigorous theory in Section [2] in a limited case, and il- 
lustrated in a relatively formal approach to an example in this section. In 
the analysis here the major issue was purely the algebraic complexity. Thus 
the main outcome of this section is to empower others to analyse any of the 
broad class of stochastic reaction-diffusion pdes identified in section [lT] that 
may be of interest in specific applications. 



4 Vary the origin of space-time to improve 
modelling 

Here reconsider linear diffusion with stochastic forcing 

Ut = u^x + 23(t,x), (39) 

on an infinite spatial domain in one dimension and for some yet to be de- 
fined noise process 23. Indeed the analysis of the next two sections almost 
always uses classic calculus and so also apply to deterministic, time depen- 
dent, forcing ?B as well as to stochastic forcing. Following a compact release 
of material it is natural to place the origin of the spatial coordinate system 
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at about the location of the release. Then the spread of material over time 
will be approximately symmetric about the spatial origin. In deterministic 
diffusion-based systems Suslov and Robert [H] showed how to choose the 
space origin for a given compact release: the optimal choice eliminates the 
slowest t~ ^-transients in the approach to the self similarity solution. 

However, when material is stochastically added/removed/moved over time, 
a marked asymmetry in the distribution about the origin generically devel- 
ops dynamically: indeed, in stochastic dynamics such asymmetry is likely 
to be a random walk, as confirmed below, that grows stochastically like 
Here we allow the reference point of the stochastic self-similarity to evolve in 
time to cater for such overall movement of the effective origin in space of the 
self-similar regime. The proposed choice removes the longest lasting memory 
integrals in the stochastic self-similar solution. 

Analogously it is expedient to change time. In deterministic diffusion- 
based systems, Suslov and Robert |44j showed how choosing the origin in 
time empowers one to eliminate the next slowest t~^/^-transients in the ap- 
proach to self similarity by better matching the variance. Correspondingly, 
for stochastic systems we seek to remove the corresponding memory integrals. 
But, instead of expressing the adaptation as a change in the origin of time, 
here we allow the relationship between real time and effective self similarity 
time to vary dynamically. 

Let's generalise the log-time transformation ^ [48j. Here scale the solu- 
tion and space-time by 

x-X(T) 1 
r = logT, e = ^^, t = t{T), u=-=uiT,0. (40) 



Then X gives the effective centre in space of the spreading material at any 
time, and T defines a pseudo-time that accounts for modifications to the 
effective width of the spreading material; we expect T ^ t to some level of 



approximation. Under the coordinate transformation (40), partial derivatives 
become 



d _ 1 d d _T' 

di ~ VTdi dt~T 



We reserve overdots for the derivative d/dr so use X' = dX/dT and T' 
dT/dt = l/t'. Substituting into the spde ^ gives 



rp> rpf 



22^3/2 2^3/2 



T3/2 c}^2 



+ 23(t,x) 
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Rearranging, using the definition ([7| of operator £, gives tlie spde in simi- 
larity variables as 

Ur = CU + X'Vfui: + {f - l)Ui:l: + t'T^/^^{t, x) . 

Now assume that the original noise process is such that t'T^^'^^{t, x) = 
W{t,^) for some cylindrical Q- Wiener process W{t,C,). Then the spde in 
similarity variables becomes 



Cu + X'Vfu^ + it' - 1)U!.^ + W. 



(41) 



Write in the Hermite basis Use the cylindrical expansion (g for the 
noise W and the spectral expansion (32) for the similarity field u{t,^). The 



properties of Hermite polynomials imply that the basis function derivatives 
Cfc^ = —ek+i/V2 and Ck^^ = 6^+2/2. Then equating coefficients of Ck, the 



linearised version of the above spde (41) becomes the component SDEs 



\kuk 



V2 



In particular, the first three SDEs are 



U2 



-U2 



V2 



(42) 
(43) 
(44) 



As in section pi solving the first SDE (42) gives the random walk of the 



amplitude uq = mo(0) + boWo{T). However, interesting results become clearer 
in this section by making bo = , the case of conservative noise, so that 
Uq = a is constant. 



Eliminate the leading transient With static X the solution of the sec- 



ond SDE (43) would involve fiuctuations generated by the memory convolu- 
tions e~'^^'^-kwi. It is these fiuctuations we remove by varying X. Let's explore 
the case Bq = and Uq = a . Then setting —^X'y/Tuo + biWi = elimi- 



nates the forcing of Ui in the second SDE (43). That is, set the spatial 'origin 



X = \p2h\ J UQ^y/T dwi{T). Because r = logT, then dwi{T) = dwi{T) / \/T . 



Thus remove the forcing of ui in (43) by choosing 



X = X(0) + V2bi [ — dwi{Ti) = X(0) + 
Jo ""o 



V2bi 



wiiT). (45) 
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Additionally choosing the initial condition X{0) appropriately [44J, then 
causes the mode Ui = for all time. With these choices for the dynami- 
cally varying spatial 'origin' X{T) of the similarity transform we eliminate 
both the transient in ui of relative magnitude 1/a/T and we eliminate all fluc- 
tuations in the Ui mode. This recognises that the noise moves the effective 
centre of the material. 



Eliminate the next transient With a fixed t = T the third SDE (44) 
would involve fluctuations e~'^'kW2 that we now remove by varying the relation 
between time and pseudo-time. Setting —-^X'y/Tui + ^{t' — l)uo + b2W2 = 



then eliminates the forcing of U2 in the third SDE (44). Assume we chose 
the 'space origin' X{T) so that Mi = . Then rearrange to dt/dT = t' = 
1 + 262^^2 /""o • Consequently, using dwi{T) = dwi{T)/\/T , remove the forcing 



of U2 in (44) by choosing real-time 



T 



t = t{0)+ I TdT + 2b2 I —dw2{T)=t{0)+T+— I ^T^dw2{T{). (46) 

Uq a 



Figure |2j shows five realisations of an example of the relationship (46 ) between 
real- and pseudo-time. Additionally choosing the initial time t(0) of the 
coordinate transform appropriately then causes M2 = for all time. 
With these choices for the 'origin' X{T) of the similarity transform, and the 
evolution of the pseudo-time T we eliminate both the transients of relative 
magnitude 1/a/T and 1/T, and all the fluctuations in the Ui and U2 modes. 

The pseudo-time is only very roughly linear in real-time (Figure [2]). Al- 
though real-time t does dominantly grow directly with pseudo-time T, the 
stochastic integral provides significant fluctuations to the relationship. The 
expectation E ^/J\ dw2(Ti) = so that we deduce E(t) = t(0) + T. How- 
ever, Ito's isometry shows 

Var/" ^idw2{Ti)= f TidTi = \T'^. 
Jo Jo 

That is, the size of the fluctuations in the pseudo-time about real time grows 
linearly in time. 

This ability to optimally move the 'origin' of space-time should also apply 
to stochastic self-similarity diffusion-like problems involving nonlinearity or 
multiplicative noise provided these effects are small enough perturbations. 
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3 4 5 6 7 

pseudo-time T 



10 



Figure 2: five realisations of tlie integral (46) for the real-time versus pseudo 



time relationship with parameters t(0) = 1 and 262/a = 0.5 . 

5 Stochastic advection and exchange in fluc- 
tuating time 

Turbulent mixing in fluids is vitally important and yet still incompletely 
understood. One route to understand such turbulent processes is to explore 
mixing in a prescribed, but random, shear flow. Majda |32] begun exploring 
a model of mixing in a flow of a linear shear multiplied by a Gaussian white 
noise process. Majda, McLaughlin, Camassa et al. [331 [HI HSI continued 
exploring aspects of the stochastic aspects of the mixing. They focussed on 
initial conditions which are statistically stationary in space. In contrast, here 
we look at the problem of mixing in a stochastic shear flow from a compact 
release in space. 

As in section |4] we find a fluctuating pseudo-time naturally arises in the 
emergent stochastic self-similarity. The 'mean' self-similarity displays clas- 
sic diffusive growth to correspond to eddy diffusivity models of turbulence. 
However, due to the stochastic shear flow, modelling fluid eddies, the concen- 
tration fields often partially reconstitute earlier times. Such reconstitution is 
an anomalous diffusion which here we resolve via a fluctuating pseudo-time 
in the stochastic self-similarity. 

But to make progress in this first treatment we simplify the shear layer 
even further than Majda, McLaughlin, Camassa et al. [331 [H], [13] . As shown 
in Figure [3| we simplify by having just two layers, or 'pipes' in which the 
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Ui{t,x) 

► ► ► ► ► > ► ► ► 

f|~|l exchange ±w(t) stochastic advection 

Figure 3: schematic diagram of a stochastic shear flow in two 'pipes' carry- 
ing some material with concentrations Uj{t, x) which exchanges between the 
'pipes' and is advected by a white noise velocity ±w{t). 



advecting velocity field is formally Gaussian white noise ±w{t), equal and op- 
posite in each 'pipe'. Notionally, these pipes correspond to near neighbouring 
streamlines in a turbulent fluid flow. Some material of concentration Uj{t, x) 
in each 'pipe' is exchanged between the 'pipes' at a non-dimensional rate 1/2 
— analogous to diffusion across streamlines in turbulent flow. The governing 
stochastic pdes for this system, in the Stratonovich interpretation, are thus 
the stochastic advection-exchange equations 

— = ,iu2-u,)-wit)—, — = ,(u,-u,)+Mt)— . (47) 



Now let's explore the stochastic self-similarity modelling of these dynamics 
following a compact release of material around x = at some initial time. 

The self similarity appears with independent variables of stretched space ^ 
and log-pseudo-time r defined by 

^ = x/\/T, r = logT, T = rjib , fj = —r] + w (48) 

(using slightly different notation in that in this section overdots denote time 
derivatives d/dt). Because of the symmetry in the stochastic advection, there 
is no need to seek a moving origin in space x that we addressed in the pre- 
vious section |4j The pseudo-time T is defined in terms of the auxiliary 
Ornstein-Uhlenbeck process ?7, that is then combined with the forcing again 
to drive T (interpret in the Stratonovich sense). Curiously, the T?7-system is a 
well-known irreproducible 'kernel' in stochastic slow manifold analysis: Chao 
and Roberts [HI §4] identified that the Try-system is its own stochastic slow 
manifold model; they analysed the corresponding Fokker-Planck equation to 
argue that on long time scales dT = \dt -\- -^dwi for effectively indepen- 
dent noise dwi. This long time scale SDE accounts for the mean growth of 
pseudo-time T with time t, and the fluctuations thereon, as shown by Fig- 
ure |4j There is a good physical reason for the 'reversals' in T(t) displayed 
by Figure |4} differential advection by a 'turbulent' fluctuation of a lump of 
material at some location spreads material; then a reversal of the fluctuation 
reconstitutes much of the earlier distribution. Reversals in T{t) represent 
such reconstitution. 
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0-1 ' 1 — ■ 1 — ■ 1 ' 1 ' 1 — ■ 1 — ■ 1 ' — I 

5 10 15 20 25 30 35 40 



real-time t 

Figure 4: five realisations of the pseudo-time T as a function of real-time t. 
The 'reversals' in the relationship correspond to reversals in advection par- 
tially reconstituting the distribution of material at an earlier time. 



In addition to transforming independent variables, we transform depen- 
dent fields to the mean and difference variables 



(mean) 
(difference) 



The mean is scaled to decay like 1/vT as in classic dispersion, but the 
difference field is scaled to decay faster, like 1/T, and moderated by the 
Ornstein-Uhlenbeck process 77. In modelling turbulent mixing, the macro- 
scopic dynamics of the mean field is of prime interest. 

Stochastic self-similarity emerges from the stochastic 'turbulent' mixing. 
Straightforward change of variable algebra derives that derivatives 



di 



TjW 



and 



dx 



]__d_ 



Substituting these transformations into the advection-exchange pde (47) and 



rearranging, the transformed system for mean and difference variables in log- 
pseudo-time is 



du 



1 A ^du 

du _^ 2 



dv 

dv , ^dv 



dr 



(49) 
(50) 
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The second spde (50) indicates that exponentially quickly in log-pseudo- 
time r, the difference field v — ?■ —du/d^. Substitute this limit into the first 



PDE (49) and it becomes 



du 



Cu 



du d'^u 



(51) 



This well known linear PDE thus describes the emergent dynamics of the 



advection-exchange system (47). 



The solutions of the PDE (51) settle on a Gaussian oc G{^). Writing the 



field m(t, ^) in the spectral expansion (32), and then equating coefficients in 
the PDE (51) leads to the system itk = —\kuk , k = 0, 1,2, ... . Hence all 
modes tend to zero like e""^/^ or quicker, except the k = mode which is 
constant. Thus u = aG{0 + 0{e-^/'^) for some constant a. The attraction 
of these equilibria in the transformed variables predicts the emergence of the 
classic self-similar spread in physical variables that the mean concentration 
l(ui + ua) = aT-^/'^G{x/VT) + 0{l/T) as t ^ oo , albeit expressed in 
fluctuating pseudo-time. 

This generic emergence of the spreading Gaussian is an appealing paral- 
lel with turbulent eddy diffusivity models. The fluctuations in the stochastic 
self-similarity emphasise the difficulty of just one part of traditional deter- 
ministic models of turbulent mixing. In real turbulence the situation is vastly 
more complex in that we picture many stochastic eddies occurring together. 
In such a situation the fluctuations that each may generate individually may 
well average to a smaller net effect. Further research is needed. 



6 Conclusions 

We demonstrated that centre manifold theory provides a straightforward and 
rigorous way of deriving the functional form of similarity solutions of nonlin- 
ear stochastic diffusion, and proving the emergence of stochastic similarity 
from quite general compact initial conditions. In particular, sections [2] and [3] 
derived explicit results for a stochastic Burgers' equation and a stochastic 
cubic reaction-diffusion equation. Section |4] then showed that we could vary 
in time the location of the origin in space, and vary the rate of time, to opti- 
mally describe the stochastic self-similarity. The last section [5] then used an 
analogous stochastic similarity transform to illustrate the emergence of an 
anomalous eddy diffusion process in a toy turbulent mixing problem. These 
techniques appear promising for useful modelling a wide class of stochastic 
nonlinear diffusion-like problems. 
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The analysis of the last two sections |4] and [5] involves purely classical 
calculus and so also applies to deterministic w{t). 
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